Finite-temperature properties of the Hubbard chain with bond-charge interaction 
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We investigate the one-dimensional Hubbard model with an additional bond-charge interaction, 
recently considered in the description of compounds that exhibit strong ID features above the 
temperature of ordered phases. The partition function of the model is exactly calculated for a value 
of the bond-charge coupling; the behavior of the specific heat and spin susceptibility as a function 
of temperature is derived at arbitrary filling, and particularly discussed across the occurring metal- 
insulator transition. The results show that the bond-charge terms weaken the spin excitations of 
the system. 

2001 PACS numbers: 71.10.Fd; 71.27.+a; 71.30.-fh; 05.30.-d 



I. INTRODUCTION 



In recent years the discovery of materials that in some 
energy regimes exhibit a strong one-dimensional charac- 
ter has renewed the investigation of models of interacting 
electrons in low dimensional lattices. Within this con- 
text, an increasing interest is nowadays devoted to the 
effects of unconventional correlation mechanisms, differ- 
ent from the usual charge-charge interaction terms be- 
tween electrons on the same site {U) and on neighboring 
sites {V). In particular, models have been considered 
10-^ which also account for the modification of the elec- 
tron hopping motion by the presence of particles with 
opposite spins (correlated hopping); such kind of terms 
are also called bond-charge interactions, since they actu- 
ally describe the interaction of charges located on bonds 
with those that are located on the lattice sites. 
The first field of application of such kind of models in 
condensed matter was the description of 7r-electrons in 
conducting polymers such as polyacetylene (CH)^; in 
particular it has been found Q that, according to the 
strength of bond-charge coupling the dimerization of the 
polymeric chain can be enhanced or destroyed. 

More recently, bond-charge models have been consid- 
ered 1^ to explain the rich temperature-pressure phase 
diagram observed for the Bechgaard salts Q, i.e. the lin- 
ear chain organic compounds such as (TMTSF)2X and 
(TMTTF)2X, where CIO4 or X = Br. Indeed, for these 
materials it has been noticed that the spin density wave 
(SDW) and the superconducting (SC) phases are adja- 
cent, the symmetry of the SC order parameter being of 
p-wave character, rather than d-wave like in cuprates. 
Recent studies have suggested that the presence of bond 
charge terms (whose coupling constants may depend on 
the pressure) could explain the interplay between the p- 
wave SC and the SDW orders, with varying pressure and 
electron density p (filling). 

This idea has yielded a remarkable effort to the investi- 
gation of electron models with correlated hopping terms, 
through both analytical and numerical methods. In par- 
ticular results have been obtained for the ground state 



as well as for the low-temperature limit, by the 
method of Bosonization ||l|,^ , which has allowed to sketch 
out the phase diagram with respect to the bond-charge 
coupling constant. 

Despite such remarkable results, a satisfactory compari- 
son with experimental data on Bechgaard salts has not 
been achieved yet. This is mainly due to the fact that, at 
very low temperatures, 2D and 3D couplings between or- 
ganic chains become relevant (see for instance the phase 
diagram in [^); this explains, in pass, the occurrence of 
ordered phases in these compounds. As a consequence, a 
one-dimensional picture for these materials is reasonable 
only above some reference temperature (of the order of 
10^ K), which of course also depends on the pressure. 
In order to compare theoretical results with experimen- 
tal observations on Bechgaard salts, it is therefore nec- 
essary to examine such models at higher temperature, 
or to test whether the low-temperature range tractable 
through Bosonization has a non-vanishing overlap with 
the ID region of the P — T phase diagram of such mate- 
rials. 

The purpose of the present letter is to investigate the 
properties for a model of bond charge interaction at finite 
and arbitrary temperature. In particular, for one value 
of the bond-charge coupling, we shall derive the exact 
behavior of thermodynamic observables such as the spe- 
cific heat and the spin susceptibility, and discuss how the 
correlated hopping terms affect the obtained shape. 



II. MODEL: SPECTRUM AND PARTITION 
FUNCTION 



The bond-charge model we shall discuss reads as fol- 
lows: 

n^-t [1-x {h,s + njg)] 4^Cj„ + (1) 

L L L 

+ 11^ n^^h^i (nj| - n^) + ^ ^ {fi^^ + n^) 

i— 1 i—1 i—1 
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In (|^) cj^, are fermionic creation and annihilation op- 
erators on a one-dimensional chain of (say) (TMTTF)2X. 
Let each site i represent an adequately chosen unit cell 
[0 and L be the total number of sites; a =T,i is the 
spin label, a denotes its opposite, n^^ — c^j^Cj^ is the 
electron charge with spin cr, and {i, j) stands for neigh- 
boring sites. At each site i four possible states are pos- 
sible, which we shall denote as follows: | t)i = c|^||0), 
|i). = ctjO), |0). = |0), m = clcl^\0). 

The term in the first line of (Q) is the hopping term, and 
in particular the parameter x represents the bond-charge 
coupling constant (for a; = the ordinary Hubbard model 
PI is recovered); the three terms in the second line re- 
spectively describe the usual on-site Coulomb repulsion, 
a possible coupling to an external magnetic field, and the 
chemical potential. Notice that, similarly to [§j^J^ we do 
not consider here the neighboring site charge-charge in- 
teraction V, since its presence can be accounted for -in a 
first approximation- through a renormalized value of U. 
The model (|l|) is rather general, and is expected to cap- 
ture the main effects of bond-charge terms. In the present 
work we shall provide exact results for the value: 

X = 1 (2) 

For this value of the coupling constant and for zero mag- 
netic field [h = 0), the exact ground phase diagram as 
a function of U and the filling p was obtained in and 
1^ for open and periodic boundary conditions respec- 
tively. This result was derived noticing that, for x — 1, a) 
the term in U commutes with the hopping term; b) the 
hopping term naturally allows a separation of the four 
possible states defined above into two groups, namely 
A = {It), li)} and B = {|0), |iT)}. In fact, for such value 
of the bond-charge coupling, the hopping term actually 
permutes A states with B states only, but not A (or B) 
states between themselves. 

In |l^] these arguments have been generalized. In the 
first instance, further commuting terms, other than the 
on-site Coulomb repulsion, such as the magnetic field, 
can be added to the hopping part. Secondly, it has been 
pointed out that the properties of the hopping term of 
(|l|) are shared by a whole subclass of Extended Hubbard 
Models, to which the Hamiltonian ([^) belongs. More 
explicitly, each model within this subclass identifies a 
specific number and set of Sutherland Species, i.e. the 
groups of states such that the Hamiltonian only permutes 
the states related to different species, leaving unaltered 
neighboring states that belong to the same species; for 
this reason the models of this subclass have been termed 
'generalized permutators' | [T^ . The number n of Suther- 
land Species is by definition not greater than the number 
of physical states (4 in the case of a single orbital). For 
model (0) the Sutherland Species are 2, and precisely A 
and B. 

Recognizing that a model identifies (up to some com- 
muting terms) a set of Sutherland Species greatly sim- 
plifies the calculation of the partition function. The cru- 



cial point which allows that is the use of open bound- 
ary conditions, instead of the customary periodic ones; 
although in the thermodynamic limit the bulk proper- 
ties are not affected by either choice, the calculations are 
more straightforward for the former. Indeed in an open 
one-dimensional chain the set of eigenvalues of a gener- 
alized permutator is equal to that of an ordinary permu- 
tator between n objects, i.e. the effective dimensionality 
of the Hilbert space is reduced (reduction theorem). As 
a consequence of that, the degeneracy of the eigenvalues 
can also be computed, simply counting the ways one can 
realize a given configuration of Sutherland Species. Such 
observations are rather general and have been used, for 
instance, to derive the exact thermodynamics of an Ex- 
tended Hubbard model of the above subclass [|3|. We 
shall apply them here to obtain the partition function of 
the bond-charge model (|l|) for the value (||). 
The reduction theorem is proved when realizing that, ac- 
cording to what observed above, the relative order of any 
sequence of states belonging to the same species is pre- 
served ; the Hamiltonian can therefore be diagonalized 
within each subspace of given set of sequences. In the 
case of model (|^) , the sequences Sa and Sb of A-species 
and -B-species are separately preserved. Moreover, since 
in this case one has two Sutherland species, each invari- 
ant subspace is in a one-to-one correspondence with the 
states of a spinless fermion (SF) space; if the A species 
relates to occupied and the B species to empty sites of 
the SF space, the form of an effective Hamiltonian for 
the SF problem is that of a tight-binding model, for the 
latter can be regarded to as a permutation between oc- 
cupied and empty sites. The eigenvalues therefore read 
— 2i ^^j^ cos /c n^, with quantum numbers valued 
or 1, and k — nn/ [L + 1), with n = 1, . . . , L. The num- 
ber of SF equals that of A-species objects (Na) and one 
has J2k "fe ~ number of empty sites is then 

Nb = L — Na- Under open boundary conditions, the 
specific sequence is irrelevant to the action of the hopping 
term, so that all the subspaces share the same spectrum; 
in general this does not hold under periodic boundary 
conditions (see for instance or jl^). The Fock space 
is thus reorganized in terms of states defined by species 
A and B and their related degeneracy. 
The inclusion of the further (commuting) terms sim- 
ply lifts the degeneracy of the eigenvalues of the hop- 
ping term, yielding a spectrum which depends, apart 
from {n^}, also on Ni-^ and on N^, the latter be- 
ing the eigenvalues of the operators X^i^i'^iT^U 
12f=i f^iti^ ~ f^il)- The spectrum of model (Q) thus reads 

E = E{{n^};N^r,N^)^ (3) 

k 

where ek = —It cos fc; the identities — = 2N^ — Na 
and N ^ + Ni = Na - 2iVx| have been exploited. In 
the case of zero magnetic field {h = 0) one recovers the 
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spectrum that was minimized in |7j at fixed number of 
particles (/x=const=0) to obtain the ground state phase 
diagram U vs p. 

The degeneracy g corresponds to the different ways one 
can realize a configuration of Sutherland species, with the 
constraint that the total numbers and N-^ appearing 
in ^ remain unchanged; a simple calculation yields 

.(E(i,.in;A'.,;A',))^ (5;;) (4) 

The rearrangement of the Fock space deriving from 
the identification of the Sutherland Species allows a 
straightforward calculation of the (gran-canonical) parti- 
tion function 



L-Na Na 



L 



(5) 



(1 + e'^^^"^))^ n + ^l-P{e^,-^.'^^,,|3,u,h))]^ 



k=l 



In the second line of (ph we have defined 



fi*{fi,P,U,h) = ^ + iln — 



2 cosh /3 ft, 



exp 2/3(^ - [//2) 



(6) 



/3 = \/{kBT) being the inverse temperature. Notice also 
that the product over k resulting in (|^) is in form simi- 
lar to the partition function of a tight binding model of 
spinless fermions, where /i* plays the role of an effective 
chemical potential renormalized by the interaction U , the 
magnetic field h and the temperature itself. 



III. RESULTS AND DISCUSSION 

By means of the partition function derived in 
the previous section, one can calculate the thermody- 
namic observables from the gran potential (per site) 
Lo — — limi^^oo fc^TlnZ. In doing that, we have elim- 
inated the chemical potential /i in favor of the filling 
through the relation p = dm /dp, as usual. 

In fig.0 we have plotted the specific heat (per site) Cv 
as a function of the temperature. In particular, in the 
top figure we have examined the case of half filling (i.e. 
p = 1) and zero magnetic field [h = 0), for different 
values of the on-site Coulomb repulsion U . One can ob- 
serve that, across the value U/t = 4, the low-temperature 
behavior of Cy changes from linear to exponential; ex- 
plicitly, for U < At we have 



Cy with 7 = 



whereas for U > At 



fc|7r 



6tv/(l-(f/W) 



, (7) 



Cv ~ ks 



{u-Aty 



(U-it) 



(8) 



This is a finite-temperature effect of a metal-insulator 
transition, in accordance with the result obtained in [Q,^, 
where a charge gap Ac ^ U — At is shown to open in the 
ground state for U > At. We recall that for x = (i.e. 
for the ordinary ID Hubbard model) no metal-insulator 
transition occurs; the bond charge term thus seems to 
give rise to a finite critical value Uc, increasing from 
to At as the coupling x is varied from to 1. It is also 
worth emphasizing that such effect is opposite to the case 
of higher dimension, where the bond-charge interaction is 
found Q to lower the critical value of the metal-insulator 
transition (Gutzwiller approximation, exact in the limit 
D oo). 
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FIG. 1. The specific heat as a function of temperature. 
Top: p = 1, h = 0: the metal-insulator transition is revealed 
through the change in low-temperature behavior from linear 
to exponential across the critical value U = At. 
Bottom: Cv for different filling values: p = 0.25 (dashed), 
p = 0.50 (dot-dashed); p = 0.75 (dotted) and p — 1 (solid); a 
low-temperature sharp peak emerges for non-vanishing mag- 
netic field. 

Notice that, the ratio U/t is expected to scale inversely 
with the pressure, since the increase of the latter roughly 
enhances the hopping amplitude; as a consequence, the 
passage from a metallic to an insulating state with in- 
creasing U/t is in accordance with the qualitative fea- 
tures of the Bechgaard salts phase diagram j|] . 
In the bottom fig.|l|, Cy is plotted for different filling val- 
ues, fixed ratio U/t = 2 and magnetic field h/t = 0.01. 
A sharp low-temperature peak, located at ksT ~ h, 
is observed to emerge as soon as the magnetic field is 
turned on. Interestingly, the peak becomes basically fill- 
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ing independent as p enters the range [p, 2 — p] , with 
p = cos^^ {—U / At) / -K . This amounts to the fact that, 
within this range of p, particles can be added to the sys- 
tem only in form of singlet pairs, in accordance with the 
features of the phase diagram in 0. 
In addition, one can show that 



lim lim Cy IT ^ lim lim Cy IT 



(9) 



differently from the ordinary Hubbard model, where the 
two limits are interchangeable . At half-filling and for 
\U + \2h\ \ < At, for instance, one has Cy ^ jT with 



fc|(31n2 2 + 7r2) 



67Tt^{l-{{U + 2\h\)/At)^) 



(10) 



Comparing eq.(Q to eq.(^, one can realize that (||) 
holds. Similarly, the exponential behavior, occurring 
when the gap is open, is different; namely, for |C/-1- |2/i| | > 
At 



Cv - k. 



{U + 2\h\-Aty 



{kBT/t)-^^^e 



(U + 2\h\-it) 



(11) 



to be compared to eq.(||). 

In fig.|| the specific heat of model (|l|) for a; = 1 is plot- 
ted aside the case a; = (i.e. the Hubbard model) for 
strong coupling, namely U = 8t. Notice that the ordi- 
nary Hubbard model has a low-temperature peak, whose 
origin (see e.g. ||l^) is due to spin degrees of freedom; 
the latter being not gapped, the low-temperature behav- 
ior of Cy is linear in spite of the fact that a charge gap 
is present at any U > 17 1. In contrast, in model (||), 
for a; = 1 the spectrum does not carry any spin quantum 
number, due to the rich symmetry of the model ||l^; spins 
act therefore as a sort of dummy variables. Although the 
value X = 1 is a particular one, it is reasonable to expect 
that, for continuity argument, the plot of Cy for inter- 
mediate values < x < 1 lies between the two curves. 
As a consequence, we can infer that the effect of spin ex- 
citations is weakened by the presence of the bond-charge 
interaction, at least in the strong coupling regime. 
In order to have a qualitative idea concerning the Bech- 
gaard salts ||^ , the temperature range of the figures com- 
patible with the ID regime of e.g. (TMTTF)2Br is 
ksT/t > 0.3 (indeed t - 0.1 eV and U - 1 eV). 




FIG. 2. The specific heat as a function of temperature for 
model (|l|) in the strong coupling regime {U = 8t), at half 
filling and zero magnetic field. The dotted line is the case 
a; = - i.e. the ordinary Hubbard model- obtained from [16], 
and the solid line the case a; = 1, obtained from our exact 
calculations. Continuity arguments suggest that the specific 
heat for arbitrary < a; < 1 lies between these two curves. 
The low-temperature peak originating from spin excitations 
is depleted by the bond-charge interaction. 



The depletion of spin excitation is also confirmed by 
the behavior of the magnetic susceptibility, defined as 
X = p%{dm/dh)\h^o- The calculation shows that x 
coincides with p^pA/ksT, where pA (the density of 
A-species along the chain) is a regular function of T, 
plotted at half filling in fig.^. One can observe that, dif- 
ferently jl^] from the ordinary Hubbard model (a; = 0), 
the susceptibility is divergent for T — + with a Curie- 
law behavior, reminiscent of a system of independent 
magnetic moments. This also holds for any filling value, 
suggesting that, as a; — > 1, the velocity of spin excita- 
tions vanishes, independently of U and of the filling; in 
particular, no spin gap is present for a; = 1 either. 
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FIG. 3. The spin susceptibility of model (jl|) with a:: = 1 
diverges like as T — > 0. The quantity x^bT is the 

integral of the spin-spin correlation function, from fluctua- 
tion-dissipation theorem. A change in its low-temperature be- 
havior is observed across the metal-insulator transition value. 

One can now compare such results with those obtained 
through an a low-energy approach. The latter is reliable 
only when the interaction couplings and the thermal fluc- 
tuations are small compared with the bandwidth w = At] 
however, for the ordinary Hubbard model, one can (a 
posteriori) extrapolate the results concerning the forma- 
tion of charge/spin gaps to the strong coupling regime 
too. 

Applying the bosonization technique to model (|^), 
one can show (see also ||l|,||) that a) the charge sec- 
tor behaves like that of the Hubbard model, so that 
J/c = at p = 1, and b) the spin sector is gapless for 
u' = U/t + 8a; cos(7rp/2) > 0, the spin excitations ex- 
hibiting a velocity Vs = vf\/^ — u' /nvp. In particular, 
a linear contribution to Cy from spin excitations is thus 
expected for any x as long as u' > 0. 
The exact results of the present work (figs.H and H, 
and the discussion above) indicate that any attempt to 
mimic the ordinary Hubbard model extrapolating the 
low-energy approach to cc — > 1 would fail with respect 
to the spin channel, since the behavior is actually quite 
different at x ~ 1. This also holds for the charge chan- 
nel, according to the results in which show that for 
X — \ the latter is gapless at half filling for \U\ < 4t. 

In conclusion, in the present paper the finite tempera- 
ture properties of the Hubbard model with bond-charge 
interaction have been exactly derived for the value (|^) 
of the bond-charge coupling. The results concerning the 
behavior of the specific heat and the magnetic suscepti- 
bility indicate that the bond-charge interaction tends to 
suppress the spin excitations of the system. 
Previous investigations about such model were concerned 



with either the ground state or the low-temperature limit. 
In contrast, our calculations are valid at any tempera- 
ture. We emphasize that for correlated quantum systems 
exact results are very rare at finite temperature, even for 
those models that have been proved to be integrable. 
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